Experimental Test of Universality of the Anderson Transition 



o : 

<N ■ 

as: 
s ■ 

<: 

(N : 



43 ■ 

or 
i ' 
^— > ■ 

c ■ 
c3 ■ 

^ : 
cr 



> 
o 
m 

q 

od 
o 



Matthias Lopez, 1 Jean-Francois Clement, 1 Pascal Szriftgiser, 1 Jean Claude Garreau, 1 and Dominique Delande 2 

Laboratoire de Physique des Lasers, Atomes et Molecules, 
Universite Lille 1 Sciences et Technologies, CNRS; F-59655 Villeneuve d'Ascq Cedex, francf] 
2 Laboratoire Kastler-Brossel, UPMC-Paris 6, ENS, CNRS; 4 Place Jussieu, F-75005 Paris, France 

(Dated: August 22, 2011) 

We experimentally test the universality of the Anderson three dimensional metal-insulator transi- 
tion. Nine sets of parameters controlling the microscopic details of this second order phase transition 
have been tested. The corresponding critical exponents are independent (within 2a) of these mi- 
croscopic details, and the average value 1.63 ± 0.05 is in very good agreement with the numerically 
predicted value, v = 1.58. 
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In the presence of a disordered potential, the classi- 
cal diffusive transport of a particle can be inhibited by 
quantum interference among the various paths where the 
particle is multiply scattered by disorder, a puzzling phe- 
nomenon known as Anderson localization [l| . The dimen- 
sionality of the system plays a major role, which can be 
understood qualitatively from the scaling theory of lo- 
calization Q. In dimension d = 3 and above, there is 
a delocalized-localized (or metal-insulator in solid-state 
physics language) transition — known as the Anderson 
transition. An energy "mobility edge" E C) which is a de- 
creasing function of disorder, separates localized motion 
at low energy from diffusive motion at high energy. On 
the localized side, the localization length £ diverges alge- 
braically £(E) oc (E c — E)~ v , with v the critical exponent 
of the transition. On the diffusive side, the diffusion con- 
stant vanishes like D(E) oc (E — E c ) s with, according to 
the scaling theory, s = (d — 2)v A key prediction of 

the scaling theory is that the critical exponents are uni- 
versal, that is, they do not depend on the microscopic 
details of the system, such as the correlation functions of 
the disorder, the dispersion relation of the particles, etc. 
Numerical experiments on simple models [4j-|6j such as 
the tight-binding Anderson model, have confirmed this 
universality, with a non-trivial value of the critical ex- 
ponent around v = 1.58 for spinless time- reversal invari- 
ant 3-dimensional (3D) systems 0- However, there is 
a huge lack of experimental results. In this letter, we 
present accurate measurements of the critical exponent 
v in 3D, and, by varying the various experimental pa- 
rameters, test whether its value is universal. 

Anderson localization is due to interference between 
long multiple scattering paths and is thus very sensitive 
to any mechanism destroying the phase coherence of the 
wavefunction, making its experimental observation and 
characterization very difficult [f| . In the context of elec- 
tronic transport in disordered samples, electron-electron 
interaction is sufficiently important to partly invalidate 
the one-body Anderson scenario, leading to a critical ex- 
ponent close to unity In a slightly different context, 
Anderson localization of acoustic [10( and electromag- 



netic 11- III waves has been experimentally observed. 
There, absorption is a limitation, and the critical expo- 
nent v ps 0.5 estimated in [l3[ does not agree with the 
numerical simulation result 1.58. 

Cold atomic gases are conceptually simple systems 
which can be exposed to well controlled disorder. For 
sufficiently cold samples, the phase coherence time of the 
wavefunction describing the center of mass motion can 
be kept sufficiently long for atomic interferometry exper- 
iments to be routinely performed. Atom-atom interac- 
tion can be reduced to a minimum, either by usin g di - 
lute cold gases, or by using a Feshbach resonance [14j. 
Direct observation of Anderson localization of atomic 
matter waves in a one- dimensional disordered [l5| or 
quasi-periodic fl6| potential has been reported, the dis- 
ordered potential being created by the effective inter- 
action with a detuned laser beam with a random spa- 
tial profile (speckle). Very recently, observations of the 
Anderson localization in a 3D atomic gases have been 
claimed However, the interpretation of both results 
involves heuristic assumptions, and a better understand- 
ing is needed before the observation of a mobility edge 
can be asserted. 

This difficulty prompted the use of a slightly differ- 
ent system, where the disordered potential is replaced by 
classical chaotic dynamics. In the kicked rotor, a paradig- 
matic system of quantum chaos, quantum mechanical in- 
terference tend to suppress the classical chaotic diffusive 
dynamics, and to induce a phenomenon originally called 
dynamical localization, later discovered to be an analog of 
ID Anderson localization in momentum space, by map- 
ping the kicked rotor onto a quasirandom ID Anderson 
model [l8j |. The experimental realization of the kicked 
rotor with laser-cooled atoms interacting with a pulsed 
standing wave allowed the first experimental observation 
of Anderson localization in ID with atomic matter waves 
[lij . In order to observe the Anderson transition, how- 
ever, an analog of the 3D Anderson model is needed [20l |. 
Here, we focus on the so-called "quasi-periodic kicked ro- 
tor" in which three incommensurate frequencies are used 
to generate the 3D character [2l( , and shown in [22| to be 
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strictly equivalent to an anisotropic 3D Anderson model, 
a fact further confirmed by a low-energy effective field 
theory 23 1. Meticulous numerical experiments Q have 
shown the universality of the critical exponent for this 
system, v — 1.58 ±0.02, in excellent agreement with val- 
ues found in the literature for the Anderson model jij. 

An experiment based on this system has characterized 
the Anderson metal- insulator transition [24|], with the 
first experimental determination of the critical exponent 
v = 1.4±0.3. The experimental setup, described in detail 
in (22| consists of laser-cooled cesium atoms interacting 



with a pulsed far-detuned standing wave (wavenumber 
fcj, = 7.4 x 10 6 in -1 and maximum one-way intensity 180 
mW). The amplitude of the kicks is modulated with two 
frequencies uj 2 and co 3 . The Hamiltonian reads: 



H = 



E. 
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JV-l 



K COS X [1 + £ COS (ii>2t) COS (o^i)] 5(t — n) 



n=0 



(1) 

where time is measured in units of the kicking period T\ , 
space in units of (2k,L) , momentum in units of 2hki^jk, 
with k = Ahk\T\jM [M is the atom mass) playing the 
role of an effective Planck constant {[x,p] = ik) and K is 
the average kick amplitude. The kicks are short enough 
(duration r = 0.8 /is) compared to the atom dynamics to 
be considered as Dirac delta functions. If ll>2, W3, n and k 
are incommensurate, this ID quasiperiodic kicked rotor 
has a 3D Anderson metal-insulator transition, display- 
ing localization in momentum space. Compared to |24| 
and {iH , the signal to noise and the stability of the ex- 
periment have been greatly improved. Atomic momen- 
tum is measured by Raman stimulated transitions. The 
previous Raman frequency generation setup used direct 
current modulation of a master laser diode to drive the 



Raman slave lasers [25|. This system has been replaced 
by a fibered phase modulator driven at 9.2 GHz. More- 
over, 3 (of 4) master diode lasers working with an exter- 
nal cavity setup have been replaced by distributed feed- 
back lasers, extending the experiment stability from a 
few hours to several days. These improvements lead to 
much better experimental signals, making the determi- 
nation of the critical exponent much more accurate and 
reliable, opening the possibility to test its universality. 

A Sisyphus-boosted magneto-optical trap prepares an 
initial thermal state of FWHM 4 x 2hkL, much narrower 
than the final (localized or diffusive) momentum distri- 
bution, and - by time-reversal symmetry - (p(t)) remains 
zero at all time t. We can directly monitor the dynamics 
rather than rely on "bulk" quantities such as the con- 
ductance, itself related to the diffusion constant. A good 
quantity characterizing the dynamics is (p 2 (t)). For prac- 
tical and historical reasons, the atomic momentum P is 
measured in units of two recoil momenta: 
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2.89 


VE 




4,0.1 -* 8,0.8 


6.67 


1.63±0.06 
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2.89 


V7 


VT7 


4,0.1 -»■ 8,0.8 


6.68 


1.57±0.08 


C 


2.89 


V5 




3,0.435 -> 10,0.435 


5.91 


1.55±0.25 


D 


2.89 


VE 


VTs 


7.5,0 -> 7.5,0.73 


e c =0.448 


1.67±0.18 


E 


2.00 


VE 




3,0.1 5.7,0.73 


4.69 


1.64±0.08 
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2.31 


VE 




4,0.1 -»■ 9,0.8 


6.07 


1.68±0.06 
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2.47 


VE 


VTs 


4,0.1 -> 9,0.8 


5.61 


1.55±0.10 


H 


3.46 


VE 


VTs 


4,0.1 -* 9,0.8 


6.86 


1.66±0.12 


1 


3.46 


VE 




4,0.1 -> 9,0.8 


7.06 


1.70±0.12 



Table I: The 9 sets of parameters used: k, LO2 and ljs control 
the microscopic details of the disorder, K controls the ampli- 
tude and e the anisotropy of the hopping amplitudes. The 
critical point K c depends on the various parameters but the 
critical exponent is universal. The weighted mean of the crit- 
ical exponent is v — 1.63 ± 0.05. The duration of the kicks is 
r = 0.8 fj.s for sets A-H, and r = 0.96 /is for set I. 



In the diffusive regime, (p 2 (t)) increases linearly with 
time; in the localized regime, it saturates at long times 
to a value proportional to the square of the localization 
length. In |22j, it was shown that, in order to characterize 
the Anderson transition, one should measure: 



A(i) = (p 2 (t))t- 2/3 . 



(3) 



2hk L 



(2) 



A(t) scales like £~ 2 / 3 in the localized regime, like i 1 / 3 in 
the diffusive regime and is constant in the critical inter- 
mediate regime. By measuring A(i), one has thus a direct 
access to the state of the system. K is a control parame- 
ter playing the role of the energy in the Anderson model. 
The critical value K c is a mobility edge separating a local- 
ized regime at low K from a diffusive regime at large K . 
The critical exponent v can be determined from the alge- 
braic divergence \ j(K c — K) v of the localization length 
near the critical point. For fixed S,W2,W3, one can also 
chose a path in the (K, e) plane which crosses the critical 
line "at right angle", i.e. faster, which improves the accu- 
racy on the critical exponent. Rather than measuring the 
full momentum distribution to obtain (p 2 (£)), we measure 
the atomic population in the zero-momentum class Tl (t). 
As shown in [26|, (p 2 {t)) is, by conservation of the num- 
ber of atoms, proportional to IIq 2 (t). The precise propor- 
tionality factor depends on the shape of the momentum 
distribution, which has been shown experimentally (2rj| 
to vary smoothly across the transition. The quantity 
measured experimentally, A C ^ P (K, t) = HQ 2 (K,t)t~ 2 / 3 , 
has the same critical point and the same critical expo- 
nent than Eq. ([3]). In the simplest experiment, parame- 
ters k, u>2, 0J3, £ are fixed and IIo(-ftr, t) recorded for several 
values of the control parameter K , from which A exp (K, t) 
is obtained. 

As one cannot pursue the experiment for arbitrary long 
times, making the distinction between the true critical 
point where A cxp (K c ,t) is a constant and a neighboring 
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Figure 1: Typical scaling function constructed from the ex- 
perimental data measured on the quasi-periodic kicked rotor 
(set E of Tab. [I]). The tip on the right side indicates the 
critical point of the Anderson transition, the diffusive (resp. 
localized) regime is the upper (resp. lower) branch. The dis- 
persion of the points is an indication of the uncertainty on 
the measured data. 



point of small positive (localized) or small negative (dif- 
fusive) slope is difficult. To circumvent this, we use the 
finite-time scaling procedure explained in 22j, which is 
a simple extension to the time domain of the finite-size 
scaling method used in solid-state physics (27[. We as- 
sume that there exists a scaling function F such that: 



A exp (K,t) = FU(K)t 



-1/3 



(4) 



where £(K) is an unknown function to be determined. 
The idea, described in detail in (H is that £(K) 
should be chosen so that F (^(X)i -1 / 3 ) is an continu- 
ous function, safe for the divergence at the critical point, 
with two distinct localized and diffusive branches. In 
order to do so, we try to collapse the various experimen- 
tal curves A exp (K,t), corresponding to each value of K, 
in a single curve, the various £,(K) - associated to the 
localization length on the localized branch - being free 
fitting parameters. A typical scaling function is shown in 
Fig. [T]and the corresponding fit parameter function £ in 
Fig. [21 Comparison with our previously published data, 
e.g. Fig. 3(c) in [24| , shows a dramatic improvement in 
the quality of the measurements. Each point A exp (K, t) 
results from independent experiments, making the sta- 
tistical errors rather easy to evaluate. We checked that 
the dispersion of points around the scaling function in 
Fig. [T]is of the order of A(lnA) 0.05, limited by the 
experimental uncertainties, the \ 2 P er degree of freedom 
of the fit being typically slightly smaller than 1. 

From finite-time measurements, it is impossible to ex- 
tract a truly diverging £(K). Several spurious phenom- 
ena [22| are responsible for the observed cut-off: The 
dominant one is the finite duration of the experiment. 
Several sources of decoherence, such as spontaneous emis- 
sion, collisions between atoms, residual effect of grav- 
ity because the laser beams are not perfectly horizontal, 
and the inhomogeneity related to the Gaussian inten- 




Figure 2: (color on line) Red points with error bars: Charac- 
teristic length £(K) constructed from the experimental data 
(set E of Tab. [I]). At low K (localized regime), it is pro- 
portional to the localization length. At large K (diffusive 
regime), it is inversely proportional to the diffusion constant. 
The "divergence" near K c — 4.69 is a signature of a second or- 
der quantum phase transition. Several phenomena limit this 
divergence, see text for discussion. They are taken into ac- 
count by a fit including a cut-off parameter, Eq. ([5]), shown as 
a solid curve, from which the critical exponent v is extracted. 



sity profile of the standing wave contribute to the cut- 
off. We have increased the phase coherence time up to 
200-300 kick periods, in agreement with theoretical cal- 
culations [22| ■ In order to reduce uncontrolled systematic 
effects, we chose to use the same time interval - from 30 to 
150 kicks - for the analysis of all experimental data. The 
finite duration is also responsible for systematic effects: 
It tends to slightly underestimate K c , but does not seem 
to shift significantly the critical exponent. Using a differ- 
ent interval, 30-120 kicks, produces v values not differing 
by more than 0.05. 

A typical £ (K) curve, Fig. [2j displays a clear "diver- 
gence" near the critical point (increase by more than one 
order of magnitude, much better than in It is itself 

fitted to extract the position of the critical point K c and 
the critical exponent i/, using the following formula: 



l/£(K) = a\K - K c 



(5) 



where /3 is a cut-off parameter taking into account the 
various limitations discussed above. As seen in Fig. [2] 
the fit is excellent. The fitting parameter v depends 
on the range of K where the fit is performed. A too 
small range produces a large uncertainty in v while the 
quality of the fit deteriorates for a too large range. To 
avoid any bias, we have fitted all data sets in the interval 
[0.8K C , l.2K c ], for which the x 2 P er degree of freedom is 
of the order of 1. The uncertainties are calculated using 
a bootstrap method starting from the raw experimental 
data and their error bars, ending to the determination of 
v through the construction of the scaling function. For 
data sets presented, no statistically significant anomaly 
has been detected. The statistical uncertainty on K c is 
very small, at most few 10~ 2 , but the finite duration of 
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Figure 3: (color online) Experimental test of the universal- 
ity of the metal-insulator transition. The critical exponent v, 
measured for 9 different sets of parameters A-I, is universal, 
i.e. independent of the microscopic details. The error bars 
indicate one standard deviation, measured using the exper- 
imental uncertainties and a bootstrap technique. The hori- 
zontal dashed line is the commonly accepted value v = 1.58, 
determined from extensive numerical calculations on various 
models. Parameter values: See Tab. [T] 

the experiment is responsible for a systematic shift to- 
wards low K (see above). 

In order to test the universality of the critical exponent, 
we have chosen a "reference" set of parameters, noted A 
in table HI which has the same parameters used in [24[ . 
We have then modified the loi and uj^ frequencies for set 
B. We have also modified the path in the (K,e) plane, 
either by changing K only (set C) or e only (set D). In 
these two situations, the crossing of the critical regime is 
slower, making the accuracy on v significantly worse. We 
have further modified the kicking period T± , which affects 
the effective Planck's constant S. Several smaller % values 
have been used in sets E, F and G. One larger value of £ 
was used in sets H and I the difference between the two 
sets being the duration of the laser pulses. The fact that 
very close values of K c and the same critical exponent 
are obtained is a strong indication that our experimen- 
tal system is well described by the model Hamiltonian, 
Eq. Q). 

The final results are given in table U and plotted 
with their error bars (one standard deviation) in Fig. [3] 
They unambiguously demonstrate the universality of the 
localized-diffusive transition in the quasi-periodic kicked 
rotor. Moreover, all numerical values are compatible 
(within 2 standard deviations) with the best numerical 
determinations of v = 1.58 ± 0.02, both for the kicked 
rotor and the Anderson model. They all markedly differ 
from the value 1 predicted by the self-consistent theory of 
localization [2{|. The later theory is an attempt to justify 
the scaling properties using a microscopic approach: It is 
qualitatively correct and gives simple physical pictures. 
For example, it has been used to successfully pre dict the 
momentum distribution at the critical point [26]. How- 
ever, it lacks a key ingredient: At criticality, the wave- 
functions display very large fluctuations which can be 
characterized by a multifractal spectrum 0M|3l|. Huge 
fluctuations are known to affect critical exponents of ther- 



modynamic phase transitions, it is thus no surprise that 
they also affect the Anderson transition. While quantum 
phase transitions are usually considered for the ground 
state of the system, it must be emphasized that the An- 
derson transition deals with excited states in the vicinity 
of the mobility edge, displaying much richer properties. 
Especially, ultra-cold atomic gases open the way to exper- 
imental studies of the interplay of disorder, interference 
and interactions. 

We thank G. Lemarie for useful discussions, and R. 
Holliger for help with the experiments. 
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